Six - Vertex Model with Domain wall boundary conditions. 
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Abstract 



We consider the six- vertex model with domain wall boundary conditions. We 
choose the inhomogeneities as solutions of the Bethe Ansatz equations. The Bethe 
Ansatz equations have many solutions, so we can consider a wide variety of inho- 
mogeneities. For certain choices of the inhomogeneities we study arrow correlation 
functions on the horizontal line going through the centre. In particular we ob- 
tain a multiple integral representation for the emptiness formation probability that 
generalizes the known formulas for XXZ antiferromagnets. 



1 Introduction 

The six- vertex model was first introduced in |IJ . It was solved exactly by E. Lieb || and 
B. Sutherland || in 1967 by means of a Bethe Ansatz for periodic boundary conditions. 
Later the six-vertex model was studied also in the presence of several other boundary 
conditions Domain wall boundary conditions were introduced in 1982 ||. These 

boundary conditions are interesting because they allow for the derivation of determinant 
representations for correlation functions ||. It was realised recently that the six- vertex 
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model in the presence of such boundary conditions is extremely helpful in the enumeration 

MM 



of alternating sign matrices 
was calculated in |T2 . 



The bulk free energy for these boundary conditions 



In this paper we show that for special choices of inhomogeneities, one can compute the 
free energy and some correlation functions of the system. This observation might be useful 
because we expect some properties of the model to be independent of the inhomogeneities 
i.e. to depend only on the anisotropy parameter. In the simplest situation, the correlation 
functions coincide with the ones for periodic boundary conditions. We will continue 
the line of research of [13| and will primarily be interested in the emptiness formation 
probability (EFP) which was first introduced in [FJJ. In the latter paper a multiple 



integral expression for the EFP was obtained for the first time for the ground state of the 
XXX antiferromagnet. Here we will consider the EFP for more general Bethe states and 
will generalize the multiple integral expression to these cases. 



2 Quantum inverse scattering 
2.1 Setup and notation 

We consider the inhomogeneous six-vertex model with domain wall boundary conditions. 
It is defined as the six-vertex model on a M x M square lattice with fixed boundary 
conditions: arrows on the horizontal (vertical) edges are outgoing (ingoing). Furthermore, 
spectral parameters Aj and are attached to line i and column k. We choose the following 
parametrization of the usual Boltzmann weights a, 6 and c, 



a(A) 
6(A) 

c(A) 



sinh(A — 77/2) 
sinh(A + 77/2) ' 

skill?? 
sinh(A + r]/2)' 



We want to make use of the formalism of the algebraic Bethe Ansatz. For that purpose 
the Boltzmann weights are collected in a matrix L, 



L(X) 



(I 








o\ 





6(A) 


c(A) 








c(A) 


6(A) 





\o 








V 



(2) 



The monodromy matrix is then defined as an ordered product of the L-operators, 

T(A; {fi k }) = L(A - Hm) ■ ■ ■ L{\ - ^ 2 )L(X - fx x ), (3) 
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where L(X — ^k) acts on the kih factor of the physical space C 2 *' , and the auxiliary space 
C 2 . For clarity, we will in the following often suppress the explicit dependence of T on 
{lik}- As an operator on the two-dimensional auxiliary space, T(A) can be written as 

TW = (4) 



C(A) D(A) 

where A, B, C, D are operators on the physical space C 2 . The trace of the monodromy 
operator, 

T(A) = A(A) + D(A) (5) 

is the usual transfer matrix corresponding to the model with periodic boundary conditions. 
Because of the parametrization ([]]), the monodromy matrix T(A) satisfies the following 
intertwining relation, 

R(X - //) • [T(A) ® T(/i)] = [T(A) <g> T(/x)] • i?(A - /x), (6) 

The R-matrix .R is defined by, 

i?(A)=PL(A + r / /2), (7) 

where P is the permutation operator on the space C 2 <g> C 2 . Equation (|6|) embodies several 
commutation relations between the operators A, B, C and D defined in (|j). 

The fixed boundary conditions imply the following formal expression for the partition 
function ||, 

Z M ({k}, W) = (II B(A i; {/i fc }) . . . B(A M ; W) IT) , (8) 

where ||) (||)) is the state with all spins up (down). In this paper we shall be interested 
only in the case where M is even and the {Aj} are chosen as {Xi}^i N+1 = {\}iL\ for 
% = 1,...,N, where N = M/2. This allows us to rewrite the partition function in a 
convenient way [13] . First define the states, 



\N) = B(Ax) . . . B(\ N ) |T) , (N\ = (T| C(A0 . . . C(A iV ). (9) 

Let R = rifcli a t, De the flip operator on the physical space that flips all arrows. We then 
find that, 

Z M ({X i },{/i k }) = (N\R\N). (10) 

Formal expressions for correlation functions can also be written concisely in this notation. 
E.g., the probability that all arrows located at the columns fci, . . . , k n and between the 
lines N and N + 1 are down is given by, 

(iV|R.7r fcl ...7r fc jA r ) nn 



where Tx k = §(1 — Averages like (|TTJ) can be calculated using the solution of the 



quantum inverse scattering problem for the operators 7i k fll6||T7 

fc-l M 

n k = H T(m + 77/2) U(ji k + rj/2) \ [ T^ + 77/2). 



;i2) 



i=i 



Z=fc+1 



From this expression it is clear that the correlation function (|TTJ) simplifies when the h% 
are nearest neighbours. 



3 Bethe Ansatz 



A marvelous aspect of formulae ( |TT| ) and ( T2 ) is that the set of inhomogeneities can be 
chosen such that the state \N) is an eigenstate of T. In fact, there are many choices 
possible that have this property. In this section we will derive explicit expression for the 



correlation function ( TTQ corresponding to such choices. We will closely follow a similar 
derivation given in [T3|] for one particular choice of the inhomogeneities, namely those 
corresponding to the groundstate of the antiferromagnetic XXZ quantum spin chain. To 
begin with, we fix the set of inhomogeneities {Aj} to be a solution of the Bethe Ansatz 
equations, 

a(^) fr b(Xj -X k + y/2) 

where a(A) = 1 and d(A) = Ylk=i K^~^k) are the eigenvalues of operators A(A) and D(A) 
respectively on the reference state ||). Note that these equations imply that Ylj=i d(Xj) = 
1. It will be useful to rewrite the Bethe Ansatz equations in their logarithmic form, 

<p{Xj) = 71 (mod 2vr), 1 < j < N, (14) 

where the function ip is defined by, 

* W = -'' n ^-'2> KA fc -A + ,,/2) 



k=l 

iln^ + if lBf- ^" + t"M )- (15) 



°(^) fbi V sinh(r] - A + X k 

If the {Aj} are a solution of the Bethe Ansatz equations, the state \N) is a common 
eigenstate of R with eigenvalue ±1, and of T(A) with eigenvalue t(X) given by, 

N N 

t(X) = a(X)l[b- 1 (X l -X + r ] /2) + d(X)l[b-\X-X l + r ] /2) 
i=i t=i 
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N 



i(A) (l + e- i ^)n^ 1 (A,-A + r ? /2). 



(16) 



i=l 



If {Aj} obey the Bethe Ansatz equations (|13D and {£•,} are a set of parameters, then the 



following holds [15, 18 



N N 

<Tin c ^)ii B ^)iT> 

i=i 



det t' 
det7 : 



where, 



dX j 



V,, 



1 



sinh(& - Xj) 

A useful formula that we will use in the following is, 

N N 



(17) 



(18) 



det V Yl sinh(£ fe - A { ) = J| sinh(A fe - Aj) sinh(& - £ k 



(19) 



k,l=l 



k,l=l 
k<l 



From ( |10|) it follows that the partition sum is given by the norm of the Bethe state \N). 
An expression for the norm of a Bethe state in terms of a determinant is given by the 
following formula which may be obtained by specializing the in ([17]) to {A;}, 



A' 



N 



(N\N) 



(rin c ^)ri B ^)it) 



i=l 



8=1 



sinh(^) 



N 



sinh(Aj — A,- + rj) 



N 

^ sinh (A,-A j ; 
V Mi 



det <p' , 



(20) 



where, 



9? 



,(dy{\) d<p(\)\ 



V 9A,- 



9a y 



(21) 



A=Ai 



The determinant formula for the norm of a Bethe wavefunction was first conjectured by 
Gaudin |?J. Due to the complicated nature of the Bethe wavefunction a proof was not 
available till the development of the quantum inverse scattering method. The first proof 
of the determinant formula of the norm of the Bethe wavefunction for the XXZ spin chain 
was given by Korepin ||. 

More work has to be done to obtain an expression for the correlation function (|TT|). 
Here we will treat the case when the ki are nearest neighbours. Using the solution for the 
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quantum inverse scattering ([12]) and the fact that \N) is an eigenstate of T(A) and R one 
finds, 

(7T fc+ i . - -7T fc+ „) = 11* Wj+W 2 ) " • ( 22 ) 

Since + 77/2) = 0, the inverse eigenvalue t~ l (nk + rj/2) takes the simple form, 

smh(Ai - Hk + rj/ 2 ) 

The action of a product of the operators D on \N) can be calculated and is given by, 

n N N+l N+2 N+n N+n 

II D (wn B ^)it) = E E ••• E G n,.,i n ({\}tr) n b ^*)it>, 

j = l k = l il=l 12 = 1 in = l k=l 

(24) 

where the function G is given by, 

n N+l 

^ 1 ,..., 4 „({A 4 }^) = n^)c(A,-A^ + r / /2) J] b -\\ H - X k + V /2). (25) 

Z=l fe=l 

fc^ii,...,ij 

We will set \N+j = l^k+j + v/ 2 to calculate the n-point correlation function (|22|). Since 
<i(/ifc) = this means that the sums in ( p4|) only run up to = N. 

From (p4]) it is seen that we need to calculate the scalar products of the type, 



Q , JX1JX . u (_m ngr g^) n; =1 go^ + y 2 ) id 

^U A iJ) l A l7 ' ' ' 5 ^V-n,/ifc+l, • • • — (iVliV) ' 

(26) 

Using (|17D and (|T9|) one may express S 1 as a ratio of determinants, 



cv r \ i r \ \ n TT smh(Ajv_ n+ j — A 

j smh(/i fc+i - 



■X 



n n -A- sinh(A^ - Ajy_ n+i ) -pr -A- sinh(A^ - + 77/2) det ^({Aj}, {n k +j}) ^ 
x 11 sinh(A, - - 7^/2) 11 11 sinh(A, - A^^- + 77) det ^'({A,}) " " ' 



The first N — n rows of the N x N matrix ip' are equal to those of the matrix if', but the 
other n rows are different. 

^ = dj, l<7<iV-n, (28) 

V>« = — rr, anhr? 1 7T --, iV-n + l<7<iV. (29) 

J smh(Aj - /Ufc+j - 77/2) smh(A i - yU fc+i + 77/2) 



Finally, we can rewrite the ratio of the determinants in (|2"7) as one determinant by in- 
verting ip', 

g£ = det W- 1 ) ■ (30) 

To proceed we have to calculate the matrix x/i'tp'* 1 . The first N — n rows of this matrix 
can be easily calculated using (|28|), 



!<*<iV-n. (31) 



In the next section we will calculate the other rows in the limit (M — > oo) for a particular 
set of solutions of the Bethe Ansatz equations. 



4 Thermodynamic limit 



In this section we will describe the thermodynamic limit M — > oo. To be able to take this 
limit we need some information on the distribution of the solutions of the Bethe Ansatz 



equations. The solutions of (|13|) fall into two classes depending on the value of rj. These 
are the socalled massive regime, where A = coshrj > 1 and the massless regime where 
|A| < 1. In this paper we will concentrate on the massless case only. Since |A| < 1 
we will use the parametrization 7 = irj and furthermore, we will restrict our attention 
to the interval n/2 > 7 > 0. We will consider the class of solutions of (|T3|) for which 



the imaginary part of each Xj is either or n/2. These are the so called 1-strings in the 



language of In the limit M — > 00, the solutions we consider thus belong to a directed 
contour C (Figure |I|) which is defined by, 



C = (—00, 00) U (00 + in/2, —00 + in/2). 



(32) 



ivr/2 



Figure 1: The contour C in the complex plane. 

Now we will derive the logarithmic form of the Bethe Ansatz equations in a more 
precise manner than was done in for (0) . For that purpose we define the function p n by, 



Pn(A) 



2 arctan (tanh A cot 117/2) , for Im A = 0, 
-2 arctan (coth A tan 717/2) , for ImA = 7r/2. 



(33) 
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For sin 717 > this function is monotonously increasing (decreasing) on the line Im A = 
(Im A = tt/2). The logarithmic version of the Bethe Ansatz equations can then be written 

as, 

= 27m,, (34) 

where the function <p is given by, 

M N 

0(A) = 5>i(A - p k ) - ]Tp 2 (A - A,). (35) 
k=i j=i 

The numbers appearing in the right hand side of ( pM| ) are integers for N odd and half 
integers for N even. For every (half) integer {n,} there are two solutions of the Bethe 
Ansatz equations, corresponding to the two different values of the imaginary part. A 
solution thus is uniquely specified by a set of integers and a corresponding set of parities, 
where the parity of a solution is defined by, 

4 

v = 1 Im A. (36) 

7T 

Given a set of (half) integers {r^} and a set of parities {fj}, a solution \j of (|34]) is called 
a particle. A solution Ah = Ah + i7r(l — t>h)/4 to the equation, 

0(A h ) = 2nm, m g {n^ or v h g {^}, (37) 

is called a hole. In the thermodynamic limit, the particles and holes have finite distribution 
densities p p and ph, defined by, 

Mp p (A)dA number of particles in [A,A + dA], (38) 
Mp h (A)dA number holes in [A, A + dA]. (39) 

Note that the 1-form dA has a direction corresponding to that of C. The density p tot of 
the total possible solutions, or vacancies, is given simply by, 

Ptot(A)=p p (A)+p h (A). (40) 

Since we are dealing with an inhomogeneous model, it will be useful to define the densities 

Ptot by, 

1 M 

Ptot(A) = — Ptot(A - p k ), (41) 

k=l 

The corresponding particle and hole densities are given by, 

Pp(A-pfc) = ^(A)ptot(A -p k ), (42) 
Ph(A-p fe ) = (l-^(A))p tot (A-/i fc ), (43) 
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where the Fermi weight i?(A) is given by, 



0(A) = (44) 



Using these densities we can take the limit M — ► oo. It follows that (|34 ) in the thermo- 
dynamic limit can be written as, 



lim —<b(X) = 7T ( -1 + 2 / 



X Ptot (X')dX' ) , ( 1+) 



where the integration is along the contour C, Figure p]. Differentiating with respect to A 
we find the formula, 



p tot (A) = K^X) - J K 2 {X - A')p P (A / )dA / , (46) 
The function K\ ot is defined by, 

M 

K^\X) = lim - £ K,{X - a*), (47) 

and is given by, 



fe=l 



1,1 sin ny 

K n {X) = ^Pn(X) = ^ ginh(A _ . n7/2) ginh(A - . n7/2) • (48) 

The thermodynamic limit of if' is found using ( f4"5| ) and the fact that 0(A) = <p(A) mod 7r, 

^ = -2vri (M^p tot (A i ) + K 2 (A, - A,)) . (49) 

Now we are in a postion to return to the calculation at the end of the previous section. 
Remember that we want to calculate the last n rows of the matrix ip'ip'^ 1 . For that 
purpose we recall from (|29"D that, 



= -2vri ( ptot(Aj - /i fe+t ) + / K 2 (Xj - X')#(X')p tot (X' - ^ +i )dA' ) , (50) 



where in the second line we have used (0) and (f46|) . From (|49|) however and the fact that 
is symmetric it follows that this is precisely equal to, 

// 1 Ptot(A/ — /ifc+i) , / r i \ 

*'— = lS P»(A,) ^ (51) 
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We thus conclude that, 

= Sij, l<i<N-n (52) 

<** = M fto ,(A,) ' (53) 

The determinant of this matrix can be written concisely as, 

1 n 

det (W" 1 ) = det£— Hp^ t (\ N - n+j ), (54) 
where the n x n matrix S 1 is given by, 

^ij = Ptoti^N— n+j ~~ P-k+i)- (55) 

Finally, using (p2[) , fl24|) , fl2~7|) and fl54|) , the emptiness formation probability can be 
written as, 

- TV iV n 

\7Tfc+l • • • TTfc+n/ = ==f 

: E • • • E H d x ^ W» II^(AiJ, 

M J_J_smh(/i fc+i - n k+m ) h=i i n =i i=i 

l<m 

(56) 

where the function H is given by, 

HY/x i r„ ia detS({A i; },{/i fc+/ }) 



l<m 

n / l-l 



!<m 

FI I II sinn ( A ^ ~ Mfc+m ~~ i^/ 2 ) I J sinh(Ai ; - /2 k+m + ij/2) J . (57) 

1=1 \m=l m=l+l / 

The last step in deriving an expression for the emptiness formation probability in the 
thermodynamic limit is to replace the sums in ([56]) by integrals using the discussion 
above. We then arrive at the following multiple integral expression, 

(* Jh • • • * O = Yt — 7 / • • • / ^ A <>' f[ * (58) 

H smh(// fe+i - Hk+m) Jc Jc 1=1 

l<m 

where the Fermi weight i?(A) is define in (f|4]). We remind the reader that the integration 
is along the directed contour C, equation Q32|), Figure [I] 
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5 Conclusion 



In this paper we have obtained a multiple integral expression for the emptiness formation 
probability (EPF) on the central horizontal line of the inhomogeneous six-vertex model 
with domain wall boundaries. We derived this expression in the thermodynamic limit 
when the inhomogeneities are chosen from a particular set of solutions of the Bethe Ansatz 
equations, namely those without a bound states but otherwise arbitrary. This result is a 
first step to obtain an expression for the EPF for general solutions of the Bethe Ansatz 
equations, i.e. also for bound state solutions. We expect that certain properties of the 
EPF are independent of the special choice of inhomogeneities. Ultimately we hope to learn 
more about such properties by studying the EPF averaged over Bethe Ansatz solutions. 
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